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We investigate the effect of dynamical fermions on the sphaleron transition 
rate at finite temperature for the Abelian Higgs model in one spatial dimen- 
sion. The fermion degrees of freedom are included through bosonization. 
Using a numerical simulation, we find that massless fermions do not change 
the rate within the measurement accuracy. Surprisingly, the exponential de- 
pendence of the sphaleron energy on the Yukawa coupling is not borne out 
by the transition rate, which shows a very weak dependence on the fermion 
mass. 



1 Introduction. 



The 1+1 dimensional Abelian Higgs model merits interest for physical prop- 
erties it shares with the electroweak theory. In particular, it has topologically 
distinct minima of the energy, corresponding to different winding numbers 
of the scalar field. Transitions between these states are possible, at zero 
temperature through quantum tunneling, and at finite temperature also by 
thermal activation. What makes these transitions physically very interesting 
is that they are accompanied by an anomalous change in the fermion number, 
if the model includes a chiral coupling of the gauge field to fermions. In the 
electroweak theory, these transitions were, in all probability, responsible for 
the erasure of the primordial baryon asymmetry. Moreover, they may have 
led to electroweak baryogenesis. Processes of these type are called sphaleron 
transitions, owing their name to sphalerons, the lowest-barrier configurations 
separating energy minima. 

Since processes violating the fermion number involve field configurations 
which are nonperturbatively far from the trivial vacuum, the problem re- 
quires a nonperturbative treatment. A useful nonperturbative framework is 
provided by the Euclidean lattice field theory, wherever processes at zero tem- 
perature or static thermal properties are concerned. However, the fermion- 
number violating processes in question occur in real time and at a finite tem- 
perature, and thus are out of reach for the Euclidean quantum theory. The 
problem simplifies considerably only in the classical approximation. Real- 
time thermal properties of the resulting classical field theory can be studied 
numerically after lattice discretization. It has been established recently that, 
under certain conditions, the classical approximation is reliable for the pro- 
cesses in question |[|. 

There is by now an extensive body of work devoted to numerical study of 
the sphaleron transition rate in the classical approximation. Over time, the 
attention has shifted from one- dimensional models 0, |3|, |, [|, || to realistic 
aproximations of the electroweak theory in 3+1 dimensions J7L |], [II]]. Nev- 
ertheless, the one-dimensional models have not yet exhausted their utility. 
In particular, they can be used to investigate the role of dynamical fermions, 
such as those present in the standard model, in the real-time processes of 
interest. While in 3+1 dimensions these degrees of freedom resist classi- 
cal treatment, such treatment is possible in one-dimensional models upon 
bosonization. This approach was first proposed by Roberge |12| who used it 
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to study static properties of sphalerons at a finite fermion density ||T3| . Here 
we apply this approach to study the real-time dynamical evolution of the 
Abelian Higgs model coupled to fermions. 

The effect of fermions on sphalerons has been investigated in a variety 
of ways. These include perturbation theory |17|, [18|], valence approximation 
[ |19|j , and expansion in the number of fermion families fTT |. These methods 



are used to determine the sphaleron (free) energy rather than the sphaleron 
transition rate, a dynamical quantity whose determination requires a real- 
time treatment at finite temperature. Recently Aarts and Smit ^Of used the 
expansion in the large number of families for a numerical study of fermions in 
a classical Bose field background. The latter method is costly numerically and 
is yet to yield a figure for the sphaleron rate. To the best of our knowledge, 
the current work is the first calculation of the sphaleron rate to date to 
account for dynamical fermions. 

The contents of the paper is as follows. In Section 2 we discuss the 
bosonized form of the model and its vacuum structure. In Section 3 we de- 
termine the variation of the sphaleron energy with the Yukawa coupling. Our 
numerical results for the sphaleron transition rate are presented in Section 
4. Section 5 contains the discussion. 



2 The model 

Our starting point is the Lagrangian density in two-dimensional space-time 

C = ^{d.-ie^A^ -\F^ + l -(D,<P){D^y - ^ - v 2 ) 2 

- V [<P(iH> + #75^) + 4>*{^ - *$75V0]> (!) 

using the standard notation for the two-component spinor fermions ip, the 
complex scalar field (ft, and the U(l) gauge field A^. Here = — 2ieA^, 
where e is the gauge coupling. We assume the fields to satisfy periodic 
boundary conditions. The scalar self-coupling and the Yukawa coupling are 
A and y, respectively. This model is to be regulated such that the gauged 
current, ip^rfip is conserved while the vector current obeys the anomaly 
equation 

d ^^ = d,.r = -^F, u . (2) 
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In its global form, the anomaly equation means that the variation of the 
baryon number B equals that of Chern-Simons number Nqs'- 



5< B -^> = l(/ d ^ + ?/H =a (3) 

Better suited for our purposes is the Bose-equivalent form of the Lagrangian: 
+ Y (<Pe- 2i ^ lu + (f,* e 2i V^) (4) 



This form is obtained by introducing a real scalar field a related to the 
fermion currents via[] 

7 / ^ \a ev ^ a 1 

'71 \ IX 



ip7v75ip = ^[d^ - e ^A t \ (5) 

\/7T JlX 



The vector current obviously obeys the anomaly equation (Q). 

Two comments are in order with regard to (|J). First, note that, in order 
to satisfy periodic boundary conditions for ([|), the field a only needs to be 
periodic modulo \fhn. Similarly to the winding number of the scalar field, 
the periodicity mismatch of a, / dxd x a/\^%7t 7 changes by an integer under 
topologically nontrivial gauge transformations. However, unlike the winding 
number of 0, the periodicity mismatch of a cannot be changed dynamically, 
because the time derivative of a does satisfy periodic boundary conditions. 
Thus, imposing periodic boundary conditions on a is a matter of a gauge 
choice. We will assume a to be periodic in space when we solve the model 
numerically. 



1 These formulae are slightly different from the standard bosonization formulae which do 
not contain the vector potential in the bosonized expressions for the currents. The reason 
for this difference is that in our model the fermions couple to the gauge field through the 
axial rather than the vector coupling. The Wilson line factor which must be included in the 
gauge invariant regularized expressions for the fermionic currents are therefore different 
here from the standard one. The Wilson line in the local limit does not reduce to one and 
is responsible for the appearance of the vector potential on the right hand side of eq.(||). 
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Secondly, note that we have included explicitly Planck's constant, obviat- 
ing the fact that the bosonization essentially links quantum theories. In order 
to determine the loop expansion parameter for the model, we re-express the 
fields in units of v. For convenience, we also express the coordinates in units 
1/vyX. We then obtain for the Lagrangian 
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+ y Ue~ 2i ^^ + 0* e 2i V™ 2 /^ , (6) 

where now 0^ = — 2igA li , g = ej\f\ and y = Y/\v 3 . It is now evident 
that the loop expansion parameter for the model is h/v 2 . This parameter 
must be small for the classical approximation to make sense. In the absence 
of fermions v 2 /h is an overall factor in front of the action in path integral. 
With fermions included, h/v 2 appears explicitly in the Lagrangian @, in 
such a way that in the range of validity of the classical approximation the 
gauge coupling of the a field is weak. 

The temporal-gauge Hamiltonian density corresponding to (|1) is 



n = ~ 



E 2 + |p|2 + n 2 + ( dx(T _ ^H A y + \ D 



+ -(\4>\ 2 -i) 2 
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+ y Ue~ 2iv ^^ + (f)*e 2iv v^7^ , (7) 

where A is the spatial component of the gauge potential, whereas E, P and 
IT are the canonical conjugate momenta of A, <ft and a, respectively. In the 
following we solve numerically the equations of motion obtained from this 
Hamiltonian. It is easy to verify that the time evolution described by H, 
preserves the local Gauss' constraints 



d x E - gJ An - 2ig{Pcj) - P*<j)*) = (8) 

While the Hamiltonian (|7|) is best suited for numerical description of real- 
time evolution, the static properties of the model are made more transparent 
by eliminating the a degree of freedom. To this end, we perform the gauge 
transformation a — >• 0, A —> A + v ^' g (y ' an d a — > a + 2v^Jfa, where a = 
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Arg0. We also solve the constraint (|8|) for II and obtain for the remaining 
degrees of freedom 

n = \( E2 + p2 P + f + %( E ' -2^) 2 ) + ^pV -2^) 2 + ip' 2 

+ 7^- 2 A 2 + \{p 2 -l) 2 + 23>pcos(a), (9) 

where P p and P a are the radial and angular conjugate momenta of the scalar 
field, respectively. We see that the fermions induce a photon mass gyjh/nv 2 . 
In the classical regime, h/v 2 <C 1, this mass is small compared to 2g the 
photon mass induced by the Higgs mechanism. This form of the Hamiltonian 
is also used in the following for generating the canonical ensemble of initial 
configurations at finite temperature. 

It is important to understand how the fermions of the original formulation 
appear in the bosonized version. First, the baryon number in the a = gauge 
is minus the Chern-Simons number: 

B = --jdxA. (10) 

Keeping this in mind, we can analyze the vacuum structure of the theory. 
In the fermion-less Abelian Higgs model, there is an exact degeneracy of 
vacua labelled by the winding number of the scalar field. In these vacua, 
related by topologically nontrivial gauge transformations, the Chern-Simons 
number is equal to the winding number. Not so in the presence of fermions: 
the minima of the energy correspond to different values of an observable 
quantity, the fermion number. Hence, these minima are no longer related by 
a gauge transformation. Moreover, these states have different energies. 

To see how this comes about, consider the minimum of the energy @ for 
a macroscopically small winding number n <C L, where L is the spatial size 
of the system. Minimizing the static part of (|9]) with respect to A, we find 

A = — 1 —F—a'. (11) 

2 # l + 4nv2p2 

In the case y = we take an ansatz of a constant p. Using (|TT|), we determine 
that the constant a' = 2 jj !l minimizes the static energy, and therefore we have 
the relation 

n = -[l + -^—]B (12) 

4:71V 2 P 2 
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where n is the winding number of the field a. We see that close to the classical 
limit minima of the energy approximately correspond to integer values of the 
fermion number. If p is constant, so is a' = 2irn/L, whereas for the constant 
magnitude of the scalar field one finds p 2 = v 2 + 0((n/L) 2 ). Hence, the 
scalar self-coupling term in the potential approaches its absolute minimum 
in the infinite- volume limit. We then obtain for the vacuum energy 

f(B )/L= fi 4(i + A_) +0 (f)\ (13) 

In the general y ^ case a' is no longer constant for a minimal-energy 
winding solution. For simplicity let us consider the limit of the gauge and 
Yukawa couplings being small compared to the scalar self-coupling: g ^ 1, 
y -C 1. Then p does not deviate significantly from 1. In the classical limit 
h/v 2 <C 1 we again have that the baryon number is equal to the winding 
number of a. Consequently one expects that the fermions will appear in the 
spectrum as the solitons of the field a. This indeed is the case. Using (|TT|), 
we obtain the static equation for a: 

v 2 

a" + 2^(4^— + 1) sin a = (14) 

This is the sine-Gordon equation, which possesses soliton solutions. The 
mass of the fermion is equal to the energy of the (one-winding) solution^ 

/ %y 

Mf = 1 1 + Anv 2 /h 

Similarly, there exist multi-soliton solutions for any n. For a vanishingly 
small fermion density n/L the n-fermion state has the energy nMf. 

At larger y the winding solutions corresponding to fermion excitations 
look somewhat differently. It is not favourable energetically to keep p spa- 
tially constant. Instead in the region where the phase a varies between — it 
and 7r, the radial field p is significantly smaller than its value in the vacuum. 

2 The mass of the fermion in the Lagrangian eq.(|j) is given by the relation Mf — Tiyv. 
This however does not contradict the bosonized result eq.(Eq) since the relation between 
the Yukawa coupling y in eq.(Q) and the coupling Y in eq.(|]) is nonlinear. The standard 
bosonization procedure leads to Y oc y 2 (taking account of the normal ordering of the 
exponential of a), which is indeed consistent with eq.(fi"5|). 
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(15) 



This suppresses the contribution of the kinetic energy of the a field to the 
energy. At very large y when the mass of a is larger than the mass of p 
this is obvious since it is the only way to keep down the energetic cost of the 
winding configuration. 

Having discussed the static properties of the model, we can now iden- 
tify the relevant dimensional scales of the problem. This identification is 
important for the numerical study presented in the following. Consider the 
relevant length scales first. These are the fermion (soliton) size of the order 



^Jh/(87iyv 2 ) and the sphaleron size, of order one (the sphaleron configura- 
tion is discussed in some detail in the next section). The system size should 
be far above, and the spatial discretization (the lattice spacing) far below 
any of these scales. Next, the relevant time scales are the inverse frequen- 
cies 0(2g)~ 1 , 0(1), and 0(2y~ 1 ^ 2 ) for the gauge, radial scalar, and angular 
scalar modes, respectively. The time integration step should be chosen well 
below any of these scales. Finally, consider the relevant energy scales. The 
sphaleron energy is of order 1. If we wish to be in the range of validity of the 
sphaleron approximation, we must insist on the inverse temperature (3 ^> 1. 
On the other hand, for the classical approximation to make sense, the tem- 
perature must be well above the Higgs (0(h/v 2 )), the photon (0(2hg/v 2 )), 
and the a-meson (0(2hy/2~y /v 2 )) masses. For light fermions, an additional 
condition relating the temperature and the system size follows from (|13|). 
Knowing that the fermion coupling to the gauge field is small in the classical 
limit h/v 2 <C 1, we can compare this energy to the exact energy for B free 
massless fermions, hirB(B + 1)/L, whose relative deviation from (|13|) is 1/B. 
The latter is small if a typical value of B is large. At a finite temperature 
1/(3 this condition is achieved if nTi(3/v 2 <C L. 

In this work, we use, in view of these considerations, the following choices 
for the various parameters: gauge coupling g = \^2~E, Yukawa coupling 
< y < 0.080, the loop expansion parameter h/v 2 = 0.05, inverse tem- 
perature 13 < (3 < 15. We performed simulations with the values of a = 0.25 
and a = 0.125 and found no measurable lattice spacing dependence of the 
transition rate. Similarly, we varied the system length between 250 and 625, 
and observed no finite size effects. 
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3 The sphaleron 



As we know, transitions can occur between the vacua of different winding 
number. Of interest for the rate of these transitions is the energy of the 
sphaleron configuration (the lowest energy barrier separating the vacua). 

For the massless fermion case [y = 0), Roberge has shown that the 
sphaleron configuration is identical to the one in the no-fermion case if the 
fermion density n is macroscopically small [13]. 

For the massive case y ^ the sphaleron solution is not known analyt- 
ically. Instead, we determined the sphaleron energy, applying the extrem- 
ization method 1113, [T5I, ITS]. The idea is to look for static solutions of the 



classical equations of motion (including sphaleron configurations) by mini- 
mizing the sum of the squares of the right hand sides of the equations of 
motion for the fields: 



v ^{wJ ' (16) 

The sum extends over all the degrees of freedom fi where V is the potential 
part of 7i. The minimization is done using a simple relaxation procedure, 
i.e., integrating relaxation equations d T ifi = —d^V. In the following, we use 
the same procedure, with V replaced by V, to measure the integer part of 
Chern-Simons number. 

We used the method as described to determine the height of the sphaleron 
barrier separating the minima with winding numbers and 1, relative to the 
absolute minimum of V. 

The results are displayed in Figure |l|. 

As is clear from the figure, the sphaleron energy is well approximated by 
the relation 

E sph (Y) = E + bVy + cy (17) 

where a least squares fit yields E Q = 0.942762 ± 0.000015, b = 0.11919 ± 
0.00056, c = 3.744 ± 0.005. As expected, E coincides with the sphaleron 
energy for massless fermions, which, in turn, is equal to the sphaleron energy 
in the absence of fermions. 

Figure |] displays the values of the real and imaginary parts of the scalar 
field {u and v) for two configurations, for y = and y = 0.001. What 
is clearly visible is that for y = the fields passes through zero, and then 
winds slowly back, through the minimum u 2 + v 2 = 1, to the other side of 
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Figure 1: The energy of the sphaleron configuration relative to the vacuum 
energy. The pluses are the values obtained by the extremization method. 
The solid curve denotes the least-squares fit to the form a + by/y + dy. 

the Mexican hat. For y = 0.001, the field passes with equal speed over the 
Mexican hat, but immediately bends toward the minimum at u — — 1, v — 0. 

The dependence displayed in (|T7| ) can be given a qualitative explanation. 
The sphaleron energy measures the difference between the energy of the 
sphaleron configuration and the vacuum energy. The latter corresponds to 
the absolute minimum of the potential, which occurs, for small positive y, 
for a field value close to —v. For y = 0, the sphaleron configuration reduces 
to the fermionless case 



Corrections to the sphaleron energy at a finite y can be estimated as 
follows. First, the phase factor exp(i7rx/L) is affected strongly for any finite 
value of y. This comes about because, the scalar field must revert to its 
vacuum value in a finite interval, centered at x = to avoid an extensive 
(oc L) contribution to the sphaleron energy. Much like in the soliton solution 
of (PI), the size of this interval is 0{fi/STiyv 2 ) 1 ^ 2 , in the classical regime 




(18) 
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Figure 2: The real and imaginary parts of the scalar field (u and v) in the 
sphaleron configuration, for y = and for y = 0.001. 

where h/v 2 1. Likewise, the corresponding correction to the sphaleron 
energy is of the same order as My given by flT5|). As y grows and becomes 
0(h/8nv 2 ), there occurs a crossover from 0(\/y) to 0(y) behavior of the 
sphaleron energy correction, given, in this case, by the spatial integral of the 
Yukawa term in the sphaleron background. 

4 The sphaleron transition rate 

In the present work, we will be interested in temperatures that are large 
enough so that quantum effects can be ignored, so that the dynamics is 
essentially described by classical field theory Moreover, we will restrict our 
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attention in this work to temperatures that are small relative to the sphaleron 
energy, so that the transition rate is suppressed by the factor e~^ Eaph . Our 
choice of the gauge and Yukawa couplings, temperature, linear size, and 
discretization in space and time is discussed at the end of Section 2. 

We use the familiar technique to simulate sphaleron transitions || . Namely, 
we use a combination of Metropolis and heat-bath Monte-Carlo algorithms to 
draw initial conditions for real-time evolution from the canonical ensemble 
at temperature 1/(3, corresponding to the discretized version of @. Hav- 
ing generated an initial condition, we switch to the Hamiltonian (0), better 
suited for real-time evolution. In doing so, we set initially o = 0, whereas IT 
is determined by the Gauss' law. 

However, the procedure for measuring the sphaleron rate requires modifi- 
cation in the presence of fermions. The rate can no longer be defined, as it was 
without fermions, as the diffusion constant per unit volume of Chern-Simons 
number. Indeed, the energy is no longer a periodic function of Chern-Simons 
number, hence the average squared topological charge can no longer grow 
linearly with time. We use an alternative measurement method, in which 
the transitions are counted directly in a real-time simulation. In order to 
easily identify sphaleron transitions, the field configuration is subjected to 
relaxation (cooling). A configuration can be thought to represent thermal 
fluctuations in the vicinity of an energy minimum (a vacuum). The cooling 
eliminates thermal fluctuations. This is done by solving, for every field tpi, 
the relaxation equation 

d T <pi = -d^V. (19) 

The algorithm is essentially the same as the one employed to determine the 
sphaleron configurations, but now using the regular static potential V instead 
of V given by (]T6|). The resulting cooled configuration has an approximately 
integer Chern-Simons number, and transitions can be easily counted, as il- 
lustrated in Figure |^. 

We performed three series of simulations. In the first series we mea- 
sured the sphaleron transition rate in the Abelian Higgs model without the 
fermions, where the rate is already known from earlier work [J|, [|. The goal 
here was to compare the rate obtained by counting transitions to the one 
found as the diffusion constant of Chern-Simons number. The reliability of 
the former method depends on the frequency of transitions. If the system size 
is too large or the temperature is too high, the counting method is not accu- 
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rate because individual transitions cannot be resolved. We verified that the 
agreement between the two methods was excellent for all the combinations 
of sizes and temperatures reported here. 




50 100 150 200 250 300 350 400 



Figure 3: A sample time history of the Chern-Simons number, measured 
before and after cooling. 

In the second series of measurements we determined the transition rate for 
massless (y = 0) fermions. Table [I| summarizes the temperature dependence 
of the rate without and with massless fermions. 



p 


r (no fermions) 
(16.2 ±0.9) x 1(T 5 


r (y 


= 0) 


13 


(16.1 ± 1) x 


10" 5 


14 


(7.06 ±0.4) x 10" 5 


(6.68 ±0.34) x 


10" 5 


15 


(3.06 ±0.2) x 10" 5 


(3.05 ±0.21) x 


10~ 5 



Table 1: A comparative table of transition rate without fermions and with 
fermions at zero Yukawa coupling. 

As is clearly visible, there is no measurable effect of the massless fermions 
on the rate within the error bars. 
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In our final series of simulations we investigated the transition rate in the 
presence of massive fermions. Specifically, the sphaleron transition rate was 
measured using the cooling/counting method at a fixed values of /3, namely 
for (3 = 13, 14 and 15. We considered values of y between and 0.080. 

0.001 n , , , , , , , n 




1e-06 L 1 1 1 1 1 1 1 d 

0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 

Y 

Figure 4: The sphaleron transition rate as a function of Yukawa coupling at 
P = 13. The solid curve is R u , the right-hand side of (pPf) . The dashed curve 
is R as given by (pTj) . 

The results are plotted in Figures f|, [| and [| As can be seen from the 
figures, the transition rate, approximately constant as a function of y for 
y < 0.02, drops off for y > 0.02, the rate at y = 0.080 being a factor 2 
{(3 = 13) to 3 (j3 = 15) smaller than at y = 0. A dropoff of the rate as a 
function of y is exactly what we expect, as the sphaleron energy increases 
with y according to fllTf ). However, this rate dependence on y is much 
weaker than predicted by the sphaleron approximation, which would predict 
an exponential dropoff (solid line). We discuss this result in more detail in 
the following section. 
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Y 

Figure 5: The sphaleron transition rate as a function of Yukawa coupling at 
(3 = 14. The solid curve is R u , the right-hand side of (^). The dashed curve 
is R as given by (j21~l). 

5 Discussion 

We investigated the effect of dynamical fermions on the sphaleron transition 
rate, combining bosonization and the classical approximation. As discussed 
in Section 2, the resulting theory retains important qualitative features ex- 
pected from the presence of dynamical fermions. This theory is also inter- 
esting in its own right, since it helps elucidate the role of additional (other 
than the gauge field and the Higgs scalar) degrees of freedom in sphaleron 
transitions. Let us then summarize the lessons learned. 

First of all, comparing the rate in the presence of massless fermions to the 
rate without fermions, we see that, within the error there is essentially effect 
of including the fermions is a factor that is independent of the temperature 
for the range of temperatures considered (12 < (3 < 16). In other words, 
the inclusion of fermions does not alter the functional dependence on the 
temperature. Such behavior agrees very well with what one would expect in 
the sphaleron approximation for the rate, which, in the absence of fermions, 



14 



0.0001 




1e-07 







0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 



Y 



Figure 6: The sphaleron transition rate as a function of Yukawa coupling at 
(3 = 15. The solid curve is R u , the right-hand side of (^). The dashed curve 
is R as given by (plj) . 



where k is a numerical constant. In this formula the exponential is the Boltz- 
mann factor of the sphaleron configuration. As we saw in Section 3, massless 
fermions do not change the sphaleron energy. The temperature dependence 
of the non-exponential pre-factor results from the existence of zero modes in 
the sphaleron background which are not present in the background of a vac- 
uum configuration. Massless fermions do give rise to new zero-modes, absent 
in the no-fermion case. In our model, the new zero-mode is the a-meson. 
However, this mode is not associated with the sphaleron only, it exists both 
in the sphaleron and the vacuum backgrounds. Hence it cannot change the 
temperature dependence of the pre-factor. 

More interesting is the dependence of the rate on the Yukawa coupling. 
As one can see from figures |], |5| and |6], the rate is approximately constant 
for < y < 0.02, edging down slowly in the range 0.02 < y < 0.08. 



gives 




(20) 
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We are not aware of, nor do we attempt here, a rigorous determination of 
the transition rate in the sphaleron approximation in presence of dynamical 
fermions. However, in a first attempt to understand its dependence on the 
Yukawa coupling, one might expect that a reasonable first approximation to 
the rate would be one that based on the rate formula (p0|) for the fermionless 
case. The solid curve (R u ) in figures [|, [| and || denotes the right-hand side 
of 5 taking for E sph the y dependent value given by (|17D . It is obvious 
that the exponential decay as a function of y thus predicted is inconsistent 
with the numerical data. 

While we will not resolve this discrepancy in the present work, we would 
like to point out two possible effects of the fermions on the rate. 

The first one is related to the fact that the presence of fermions lifts 
the degeneracy of the vacuum as a function of Chern-Simons number to the 
energy of the created (anti)fermion, as is indicated schematically in Figure 
|5|. In the sphaleron approximation the rate is controlled by the sphaleron 
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Figure 7: Schematic representation of the effective potential as a function of 
Chern-Simons number. 



energy, or rather, the difference between the sphaleron energy and the vac- 
uum energy. This would be applicable for the rate from the Ncs = to 
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the iVcs ~ ±1 state. However, the corresponding Boltzmann factor for the 
reverse transition is different (larger), due to the energy difference Mf (cf 
QI5D), by a factor e@ M f . The same conclusion is reached by applying detailed 
balance at temperature T. 

To compare the sphaleron approximation with the measured rate, we 
should apply an appropriate averaging over the various types of transitions, 
that include creation as well as annihilation of (anti)fermions. To this effect, 
let us consider the following simplified model. Our one- dimensional system 
can be divided into g intervals of the order of the soliton size. Each interval 
can be in any of three states: empty, occupied by a soliton (fermion), or 
occupied by an antisoliton (but not both), each with energy Mf. Transi- 
tions can take place from the empty state to the soliton or antisoliton state 
with probability per unit time R u , while the reverse probability is Rd, by 
detailed balance equal to R u e" f . Direct transitions between soliton and 
antisoliton state are not possible. The total transition rate will be equal to 
the g times the (Boltzmann-)averaged transition rate per interval. The latter 
yields (2R u + 2R d e~P M f) / {l + 2e~P M f) = AR u /(l + 2e~ l3M f) (the denominator 
comes from normalization). The total transition rate is simply g times larger. 
One finds for the total rate 

R = l + 2e-^f Ru - (21) 



As Mf is proportional to y/Y, we thus find a non-trivial VJ^-dependent 
exponential correction. It is to be expected that a similar correction is present 
in the formula for the effective sphaleron transition rate. 

The rate R is plotted for comparison with our numerical data in the 
Figures |] through |6|. It is obvious that for 3^ > 0.02 R is still in disagreement 
with our numerical results. A somewhat better agreement between the data 
and R at small values of y is possibly coincidental, since the description in 
terms of localized non-overlapping solitons only applies at large y, such that 
exp (—/3y) is a small number. We conclude that the near independence of the 
transition rate on y cannot be explained away by the statistical distribution 
of fermions. 

Another possible explanation for the very weak dependence of the rate 
on y is as follows. In the zero fermion mass case one degree of freedom, 
the phase of the scalar field, is a zero mode, and we have a corresponding 
Goldstone boson. For y > 0, the degeneracy is broken, and the Goldstone 
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boson acquires a mass. However, for very small values of y, the symmetry 
remains approximate, and all values of the phase angle remain almost equally 
occupied. In particular, we can expect a high fermion density. This means 
that the Gaussian approximation around the field minimum (0 ps —1) will 
not be a good one, and one has to apply an appropriate treatment of the 
Sine-Gordon model instead. We intend to address this point in the future. 
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